N-Heterocyclic Carbene-Catalyzed Facile Synthesis of Phthalidyl Sulfonohydrazones: Density Functional Theory Mechanistic Insights and Docking Interactions.

N-heterocyclic carbene catalysis reaction protocol is disclosed for the synthesis of phthalidyl sulfonohydrazones. A broad range of N-tosyl hydrazones react effectively with phthalaldehyde derivatives under open-air conditions, enabling the formation of a new C-N bond via an oxidative path. The reaction proceeds under mild reaction conditions with broad substrate scopes, wide functional group tolerance, and good to excellent yields. The mechanistic pathway is studied successfully using control experiments, competitive reactions, ESI-MS spectral analyses of the reaction mixture, and computational study by density functional theory. The potential use of one of the phthalidyl sulfonohydrazone derivatives as the inhibitor of ß-ketoacyl acyl carrier protein synthase I of Escherichia coli is investigated using molecular docking.

The Stability of Hydrogen-Bonded Ion-Pair Complex Unexpectedly Increases with Increasing Solvent Polarity.

The generally observed decrease of the electrostatic energy in the complex with increasing solvent polarity has led to the assumption that the stability of the complexes with ion-pair hydrogen bonds decreases with increasing solvent polarity. Besides, the smaller solvent-accessible surface area (SASA) of the complex in comparison with the isolated subsystems results in a smaller solvation energy of the latter, leading to a destabilization of the complex in the solvent compared to the gas phase. In our study, which combines Nuclear Magnetic Resonance, Infrared Spectroscopy experiments, quantum chemical calculations, and molecular dynamics (MD) simulations, we question the general validity of this statement. We demonstrate that the binding free energy of the ion-pair hydrogen-bonded complex between 2-fluoropropionic acid and n-butylamine (CH3CHFCOO-…NH3But+) increases with increased solvent polarity. This phenomenon is rationalized by a substantial charge transfer between the subsystems that constitute the ion-pair hydrogen-bonded complex. This unexpected finding introduces a new perspective to our understanding of solvation dynamics, emphasizing the interplay between solvent polarity and molecular stability within hydrogen-bonded systems.

Impact of dielectric constant of solvent on the formation of transition metal-ammine complexes.

The DFT-level computational investigations into Gibbs free energies (ΔG) demonstrate that as the dielectric constant of the solvent increases, the stabilities of [M(NH3 )n ]2+/3+ (n = 4, 6; M = selected 3d transition metals) complexes decrease. However, there is no observed correlation between the stability of the complex and the solvent donor number. Analysis of the charge transfer and Wiberg bond indices indicates a dative-bond character in all the complexes. The solvent effect assessed through solvation energy is determined by the change in the solvent accessible surface area (SASA) and the change in the charge distribution that occurs during complex formation. It has been observed that the SASA and charge transfer are different in the different coordination numbers, resulting in a variation in the solvent effect on complex stability in different solvents. This ultimately leads to a change between the relative stability of complexes with different coordination numbers while increasing the solvent polarity for a few complexes. Moreover, the findings indicate a direct relationship between ΔΔG (∆Gsolvent -∆Ggas ) and ΔEsolv , which enables the computation of ΔG for the compounds in a particular solvent using only ΔGgas and ΔEsolv . This approach is less computationally expensive.

Coexistence of Electrochromism and Bipolar Nonvolatile Memory in a Single Viologen.

Viologens are fascinating redox-active organic compounds that have been widely explored in electrochromic devices (ECDs). However, the combination of electrochromic and resistive random-access memory in a single viologen remains unexplored. We report the coexistence of bistate electrochromic and single-resistor (1R) memory functions in a novel viologen. A high-performance electrochromic function is achieved by combining viologen (BzV2+2PF6) with polythiophene (P3HT), enabling a "push-pull" electronic effect due to the efficient intermolecular charge transfer in response to an applied bias. The ECDs show high coloration efficiency (ca. 1150 ± 10 cm2 C-1), subsecond switching time, good cycle stability (>103 switching cycles), and low-bias operation (±1.5 V). The ECDs require low power for switching the color states (55 µW cm-2 for magenta and 141 µW cm-2 for blue color). The random-access memory devices (p+2-Si/BzV2+2PF6/Al) exhibit distinct low and high resistive states with an ON/OFF ratio of ∼103, bipolar and nonvolatile characteristics that manifest good performances, and "Write"-"Read"-"Erase" (WRE) functions. The charge conduction mechanism of the RRAM device is elucidated by the Poole-Frenkel model where SET and RESET states arise at a low transition voltage (VT = ±1.7 V). Device statistics and performance parameters for both electrochromic and memory devices are compared with the literature data. Our findings on electrochromism and nonvolatile memory originated in the same viologen could boost the development of multifunctional, smart, wearable, flexible, and low-cost optoelectronic devices.

Trends in the stability of covalent dative bonds with variable solvent polarity depend on the charge transfer in the Lewis electron-pair system.

In general, the stability of neutral complexes with dative bonds increases as the polarity of the solvent increases. This is based on the fact that the dipole moment of the complex increases as the charge transferred from the donor to the acceptor increases. As a result, the solvation energy of the complex becomes greater than that of subsystems, causing an increase in the stabilization energy with increasing solvent polarity. Our research confirms this assumption, but only when the charge transfer is sufficiently large. If it is below a certain threshold, the increase in the complex's dipole moment is insufficient to result in a higher solvation energy than subsystems. Thus, the magnitude of the charge transfer in the Lewis electron-pair system determines the stability trends of dative bonds with varying solvent polarity. We used molecular dynamics (MD) simulations based on an explicit solvent model, which is considered more reliable, to verify the results obtained with a continuous solvent model.

The Impact of the Solvent Dielectric Constant on A←NH3 Dative Bond Depends on the Nature of the Lewis Electron-Pair Systems.

The present work aims to determine to what extent the value of the dielectric constant of the solvent can influence the dative bond in Lewis electron pair bonding systems. For this purpose, two different systems, namely H3 B←NH3 and {Zn←(NH3 )}2+ , were studied in selected solvents with significantly different dielectric constants. Based on the results from state-of-the-art computational methods using DFT, constrained DFT, energy decomposition analyses, solvent accessible surface area, and charge transfer calculations, we found that the stability of the neutral H3 B←NH3 system increases with increasing solvent polarity. In contrast, the opposite trend is observed for the positively charged {Zn←(NH3 )}2+ . The observed changes are attributed to different charge redistributions in neutral and charged complexes, which are reflected by a different response to the solvent and are quantified by changes in solvation energies.

The stability of covalent dative bond significantly increases with increasing solvent polarity.

It is generally expected that a solvent has only marginal effect on the stability of a covalent bond. In this work, we present a combined computational and experimental study showing a surprising stabilization of the covalent/dative bond in Me3NBH3 complex with increasing solvent polarity. The results show that for a given complex, its stability correlates with the strength of the bond. Notably, the trends in calculated changes of binding (free) energies, observed with increasing solvent polarity, match the differences in the solvation energies (ΔEsolv) of the complex and isolated fragments. Furthermore, the studies performed on the set of the dative complexes, with different atoms involved in the bond, show a linear correlation between the changes of binding free energies and ΔEsolv. The observed data indicate that the ionic part of the combined ionic-covalent character of the bond is responsible for the stabilizing effects of solvents.

P-Doped graphene-C60 nanocomposite: a donor-acceptor complex with a P-C dative bond.

Phosphorous-doped graphene can form a covalent dative bond with the electron acceptor, C60 molecule. On the other hand, C60 on graphene and N-doped graphene surfaces can only form vdW complexes. State-of-the-art quantum-chemical techniques have been used to characterise the nature of the P-C dative bond. A considerable amount of charge transfer from the P-Gr surface to C60 has been observed. This complex formation may enable enhancement in the optical limiting response with potential application in energy harvesting. The stability of the P-C dative bond has been assessed using DFT-D molecular dynamics simulations at 300 K for 10 ps.

Addition Reaction between Piperidine and C60 to Form 1,4-Disubstituted C60 Proceeds through van der Waals and Dative Bond Complexes: Theoretical and Experimental Study.

A combined computational and experimental study reveals the character of the C60 complexes with piperidine formed under different reaction conditions. The IR and NMR experiments detect the dative bond complex, which according to NMR, is stable in the oxygen-free environment and transforms to the adduct complex in the presence of O2. Computational studies on the character of reaction channels rationalize the experimental observations. They show that the piperidine dimer rather than a single piperidine molecule is required for the complex formation. The calculations reveal significant differences in the dative bond and adduct complexes' character, suggesting a considerable versatility in their electronic properties modulated by the environment. This capability offers new application potential in several fields, such as in energy storage devices.

Cyclo[n]carbons Form Strong N → C Dative/Covalent Bonds with Piperidine.

The newly synthesized C18 ring is demonstrated as the smallest all-carbon acceptor that exhibits strong electron acceptance. This study provides a quantum-chemical investigation of the electron-acceptance behavior of monocyclic carbon rings with a particular emphasis on C18 through the formation of a dative bond with piperidine. The results show that Cn rings form strong dative bonds with piperidine, whereas the respective van der Waals (vdW) complexes are higher in energy. The main driving force is the release of angle strain of cyclo[n]carbons caused by the change in hybridization from sp to sp2 associated with the formation of the dative bond. On the contrary, other sp allotropes, diynes, favorably form vdW complexes. Molecular dynamics (MD) simulations support the stability of the dative bond throughout a simulation of 20 ps. This opens up the possibility of stabilizing highly reactive C18 through dative/covalent functionalization.

Tuning the P-C dative/covalent bond formation in R3P-C60 complexes by changing the R group.

The P-C dative/covalent bonds formed in R3P-C60 complexes (R = OCH3, N(CH3)2, NC4H8) have been affected by the nature of the R group. The highest stabilisation (18.7 kcal mol-1) has been found in the last system. The contribution of dispersion energies in the stabilisation also varies depending on the R group. The nature of the P→C bond has been characterised using state-of-the-art quantum-chemical techniques including NBO, AIM and ELF. The P→C dative bond is significantly different from the prototype dative bonds appearing in H3N→BH3 as well as in the fullerene - secondary-amine complexes previously studied by us. The findings obtained through electron structure theory have been supported by 10 ps DFT-D MD simulations.

Structure-directed formation of the dative/covalent bonds in complexes with C70piperidine.

The combined experimental-computational study has been performed to investigate the complexes formed between C70 carbon allotrope and piperidine. The results of FT-IR, H-NMR, and C-NMR measurements, together with the calculations based on the DFT approach and molecular dynamics simulations, prove the existence of dative/covalent bonding in C70piperidine complexes. The dative bond forms not only at the region of five- and six-membered rings, observed previously with C60, but also at the region formed of six-membered rings. The structure, i.e., nonplanarity, explains the observed dative bond formation. New findings on the character of interaction of secondary amines with C70 bring new aspects for the rational design of modified fullerenes and their applications in electrocatalysis, spintronics, and energy storage.

Spin modification of iron(ii) complexes via covalent (dative) and dispersion guided non-covalent bonding with N-heterocyclic carbenes: DFT, DLPNO-CCSD(T) and MCSCF studies.

This paper reports a systematic quantum chemical study of the complexation of iron(ii) phthalocyanine (FePc) and Fe(ii)-porphyrazine (FePz) complexes with N-heterocyclic carbenes. The complex can be stabilised either in a singlet or in a triplet state depending on the choice of substituents in the same carbene core. Here, the spin crossover is invoked either via direct bonding or via a dispersion guided non-covalent interaction. Stabilization of intermediate spin of Fe(ii) phthalocyanine/Fe(ii)-porphyrazine type systems via complexation with carbenes is hitherto unknown. For calculations we have considered a wide range of density functional theory (DFT) levels, including GGA, hybrid and double hybrid functionals. The results obtained using DFT are further supported by domain-based local pair natural orbital-CCSD(T) (DLPNO-CCSD(T)) and multi-configurational self-consistent field (MCSCF) results.

An Isolated Molecule of Iron(II) Phthalocyanin Exhibits Quintet Ground-State: A Nexus between Theory and Experiment.

Iron(II) phthalocyanine (FePc) is an important member of the phthalocyanines family with potential applications in the fields of electrocatalysis, magnetic switching, electrochemical sensing, and phototheranostics. Despite the importance of electronic properties of FePc in these applications, a reliable determination of its ground-state is still challenging. Here we present combined state of the art computational methods and experimental approaches, that is, Mössbauer spectroscopy and Superconducting Quantum Interference Device (SQUID) magnetic measurements to identify the ground state of FePc. While the nature of the ground state obtained with density functional theory (DFT) depends on the functional, giving mostly the triplet state, multi-reference complete active space second-order perturbation theory (CASPT2) and density matrix renormalization group (DMRG) methods assign quintet as the FePc ground-state in gas-phase. This has been confirmed by the hyperfine parameters obtained from 57 Fe Mössbauer spectroscopy performed in frozen monochlorobenzene. The use of monochlorobenzene guarantees an isolated nature of the FePc as indicated by a zero Weiss temperature. The results open doors for exploring the ground state of other metal porphyrin molecules and their controlled spin transitions via external stimuli.

Non-covalent control of spin-state in metal-organic complex by positioning on N-doped graphene.

Nitrogen doping of graphene significantly affects its chemical properties, which is particularly important in molecular sensing and electrocatalysis applications. However, detailed insight into interaction between N-dopant and molecules at the atomic scale is currently lacking. Here we demonstrate control over the spin state of a single iron(II) phthalocyanine molecule by its positioning on N-doped graphene. The spin transition was driven by weak intermixing between orbitals with z-component of N-dopant (pz of N-dopant) and molecule (dxz, dyz, dz2) with subsequent reordering of the Fe d-orbitals. The transition was accompanied by an electron density redistribution within the molecule, sensed by atomic force microscopy with CO-functionalized tip. This demonstrates the unique capability of the high-resolution imaging technique to discriminate between different spin states of single molecules. Moreover, we present a method for triggering spin state transitions and tuning the electronic properties of molecules through weak non-covalent interaction with suitably functionalized graphene.

The X40×10 Halogen Bonding Benchmark Revisited: Surprising Importance of (n-1)d Subvalence Correlation.

We have re-evaluated the X40×10 benchmark for halogen bonding using conventional and explicitly correlated coupled cluster methods. For the aromatic dimers at small separation, improved CCSD(T)-MP2 "high-level corrections" (HLCs) cause substantial reductions in the dissociation energy. For the bromine and iodine species, (n-1)d subvalence correlation increases dissociation energies and turns out to be more important for noncovalent interactions than is generally realized; (n-1)sp subvalence correlation is much less important. The (n-1)d subvalence term is dominated by core-valence correlation; with the smaller cc-pVDZ-F12-PP and cc-pVTZ-F12-PP basis sets, basis set convergence for the core-core contribution becomes sufficiently erratic that it may compromise results overall. The two factors conspire to generate discrepancies of up to 0.9 kcal/mol (0.16 kcal/mol RMS) between the original X40×10 data and the present revision.

Conventional and Explicitly Correlated ab Initio Benchmark Study on Water Clusters: Revision of the BEGDB and WATER27 Data Sets.

Benchmark ab initio energies for BEGDB and WATER27 data sets have been re-examined at the MP2 and CCSD(T) levels with both conventional and explicitly correlated (F12) approaches. The basis set convergence of both conventional and explicitly correlated methods has been investigated in detail, both with and without counterpoise corrections. For the MP2 and CCSD-MP2 contributions, rapid basis set convergence observed with explicitly correlated methods is compared to conventional methods. However, conventional, orbital-based calculations are preferred for the calculation of the (T) term, since it does not benefit from F12. CCSD(F12*) converges somewhat faster with the basis set than CCSD-F12b for the CCSD-MP2 term. The performance of various DFT methods is also evaluated for the BEGDB data set, and results show that Head-Gordon's ωB97X-V and ωB97M-V functionals outperform all other DFT functionals. Counterpoise-corrected DSD-PBEP86 and raw DSD-PBEPBE-NL also perform well and are close to MP2 results. In the WATER27 data set, the anionic (deprotonated) water clusters exhibit unacceptably slow basis set convergence with the regular cc-pVnZ-F12 basis sets, which have only diffuse s and p functions. To overcome this, we have constructed modified basis sets, denoted aug-cc-pVnZ-F12 or aVnZ-F12, which have been augmented with diffuse functions on the higher angular momenta. The calculated final dissociation energies of BEGDB and WATER27 data sets are available in the Supporting Information. Our best calculated dissociation energies can be reproduced through n-body expansion, provided one pushes to the basis set and electron correlation limit for the two-body term; for the three-body term, post-MP2 contributions (particularly CCSD-MP2) are important for capturing the three-body dispersion effects. Terms beyond four-body can be adequately captured at the MP2-F12 level.

Prediction of neutral noble gas insertion compounds with heavier pnictides: FNgY (Ng = Kr and Xe; Y = As, Sb and Bi).

A novel class of interesting insertion compounds obtained through the insertion of a noble gas atom into the heavier pnictides have been explored by various ab initio quantum chemical techniques. Recently, the first neutral noble gas insertion compounds, FXeY (Y = P, N), were theoretically predicted to be stable; the triplet state was found to be the most stable state, with a high triplet-singlet energy gap, by our group. In this study, we investigated another noble gas inserted compound, FNgY (Ng = Kr and Xe; Y = As, Sb and Bi), with a triplet ground state. Density functional theory (DFT), second order Møller-Plesset perturbation theory (MP2), coupled-cluster theory (CCSD(T)) and multi-reference configuration interaction (MRCI) based techniques have been utilized to investigate the structures, stabilities, harmonic vibrational frequencies, charge distributions and topological properties of these compounds. These predicted species, FNgY (Ng = Kr and Xe; Y = As, Sb and Bi) are found to be energetically stable with respect to all the probable 2-body and 3-body dissociation pathways, except for the 2-body channel leading to the global minimum products (FY + Ng). Nevertheless, the finite barrier height corresponding to the saddle points of the compounds connected to their respective global minima products indicates that these compounds are kinetically stable. The structural parameters, energetics, and charge distribution results as well as atoms-in-molecules (AIM) analysis suggest that these predicted molecules can be best represented as F(-)[(3)NgY](+). Thus, all the aforementioned computed results clearly indicate that it may be possible to experimentally prepare the most stable triplet state of FNgY molecules under cryogenic conditions through a matrix isolation technique.

What Are the Ground State Structures of C20 and C24? An Explicitly Correlated Ab Initio Approach.

A new benchmark study has been performed for six isomers of C20 and four isomers of C24 using explicitly correlated methods, together with coupled cluster theory with large basis sets and DFT with advanced functionals. The relative energy trends obtained are extremely sensitive to the methods used. Combining our best CCSD(T)-MP2 difference with our best MP2 basis set limit, the dehydrocorannulene bowl is found to be the most stable for C20, followed by the cage at about 8 kcal/mol, and the ring at about 46 kcal/mol. For C24, the D3d cage is found to be the most stable isomer, followed at only a few kilocalories per mole by dehydrocoronene, and at larger separations by then octahedral cage and the ring, respectively. This makes C24 the smallest classical fullerene. The estimated residual basis set error of the estimated CCSD(T) basis set limit is conservatively expected to be ±1 kcal/mol. In general, DFT exhibits large errors for relative energies with RMSD values in the 8-34 kcal/mol range. However, among the DFT functionals, the DSD-PBEP86-D3BJ double hybrid comes close to our best ab initio results, while the ωB97X-V range-separated hybrid is in semiquantitative agreement.